Fingerprints of spatial charge transfer in Quantum Cascade Lasers 
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We show that mid infrared transmission spectroscopy of a quantum cascade laser provides clear cut 
information on changes in charge location at different bias. Theoretical simulations of the evolution 
of the gain/absorption spectrum for the A ~ 7.4 ycm InGaAs/AllnAs/InP quantum cascade laser have 
been compared with the experimental findings. Transfer of electrons between the ground states in 
the active region and the states in the injector goes in hand with a decrease of discrete intersubband 
absorption peaks and an increase of broad high-energy absorption towards the continuum delocalised 
states above the barriers. 



I. INTRODUCTION 

The performance of Quantum Cascade Lasers (QCL) 
has improved fast since the first experimental realisation 
more than a decade ago [lj. The emission wavelength 
now ranges between 2.95 and 190 fim (6.5 - 420 meV) 
[2, Q and high output powers have been reached [J, Q • 
The key operating mechanism of most QCLs is efficient 
and rapid transport of electrons from the lower laser state 
to the upper laser state of the next period. In order to 
achieve this, the laser structure is carefully designed to 
reach resonance conditions at operating bias in order to 
optimise tunnelling rates. Knowledge of the actual loca- 
tion of electron charges in the structure which influences 
the self-consistent electric potential is crucial to design 
the resonance conditions as well as transition energies. 

Most detailed previous theoretical studies were focused 
on level occupations and transport properties [f| 0, H, Q • 
In contrast, here we focus on the evolution of the spatial 
distribution of charge as a function of current in combina- 
tion with the optical spectra. In particular we introduce 
the concept of spatially resolved gain and absorption in 
order to better visualise the movement of charge. 

Transmission spectroscopy measurements have re- 
cently been preformed to characterise the state of an op- 
erating QCL [3 

and in this manuscript we present sim- 
ulations for the same device. Preliminary results focused 
on discrete transitions within the active region [ill ] . Here 
we provide a more thourough discussion and additionally 
address absorption to the continuum and spatially re- 
solved gain. Furthermore, we discuss some discrepancies 
between experiments and simulations. 

This manuscript is organised as follows; section II 
presents our mathematical formalism. Our numerical re- 



sults are presented in section III. We begin the discussion 
by analysing the gain spectrum (sect. Ilia) and introduce 
spatially resolved gain in section Illb. We then study 
high energy absorption (sect. IIIc) and discuss the dis- 
crepancies between experiments and simulation in section 
Hid. A short summary and conclusion close the paper 
(sect. IV). 
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II. THEORY 



Our device simulations are based on the non- 
equilibrium Greens functions approach (NEGF) [l2| tak- 
ing into account the full non-diagonal self-energies as out- 
lined in Ref. [l3[. The device is described by a periodic 
repetition of one section of the structure, which we refer 
to as period. We employ a set of basis states \l/ a (z)e lk ' r , 
where k and r are vectors in the (x, y) plane perpendicu- 
lar to the growth direction z to describe each period. For 
the self-consistent transport calculations we use Wannicr 
states [HI for ^ a (z), taking into account the 14 lowest 
states per period together with the corresponding states 
of both adjacent periods. For the stationary state we 
calculate the eigenstates of the free particle Hamiltonian 
including the mean field contribution, the Wannier-Stark 
(WS) states. Here we add further states (assumed to be 
unoccupied) in order to obtain a better resolution of the 
high-frequency absorption. The quasi-energies E a and 
widths T a of the WS states are extracted from our self- 
consistent NEGF simulations, where the levelshifts and 
widths of the extra states are set in accordance to the 
other high-energy states. 

The material gain spectra g{u>) are evaluated by a 
simple approach assuming Lorentzian lifetime broaden- 
ing and restricting to diagonal dephasing, see also ap- 
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where d is the length of one period, e < is the elec- 
tron charge, eo is the dielectric constant of vacuum, e r 
the relative dielectric constant of the well material, c is 
the vacuum speed of light, E a , n a are the energy and 
sheet electron density of state a, respectively, N is the 
total number of periods considered and huj is the photon 
energy. Correlations in the scatterin g en vironment can 
reduce the absorption/gain lincwidth [15l |. However, this 
effect is neglected here for simplicity. In other words, 
the width of the individual transitions is estimated by 

r a ,/3 = r Q + Tp. 



III. NUMERICAL RESULTS AND DISCUSSION 

We consider the device of Ref. [l(|, which is a lattice 
matched Ino.53Gao.47As/Alo.4sIno.52As InP-based design 
with a four well active region emitting at A m 7.4 /im 
(168 meV). The laser ridge is L — 2.5 mm long, 55 /itm 
wide, and consists of TV = 35 periods. The nominal dop- 
ing density is 9.1 x 10 10 cm -2 per period. All calculations 
are performed for a lattice temperature of 77 K and we 
use the material parameters of Ref. [16j . The conduction 
band profile and 19 Wannier-Stark states of the central 
period for two different biases can be seen in Fig.[T] Here 
we also show the energetically resolved electron distribu- 
tion, which reads 



n(E,z) 



2 (for spin) 
2^4 



£ rp(z)Mz)(-i)G<p(k,E), 

a.p,k 

(2) 

where G < (E) is the correlation function, which is 
related to the density matrix by (a^(k)a a (k)) = 

J d£JG^ /3 (k, E), where a£(k) (a Q (k)) are the cre- 
ation (annihilation) operators [l7| and A is the sample 
area. The charge is located in the injector region around 
z = 10 nm for low bias (100 mV per period), but it is 
transfered to the injector region around z = 40 nm under 
operating conditions (230 mV per period). 

The calculated and measured spectra are displayed in 
Fig. [21 The major features are four absorption peaks at 
approximately 200, 280, 370, and 420 meV as well as a 
gain peak at 150 meV at high currents. The lower en- 
ergy absorption peaks disappear with increasing current, 
while the broad high-energy absorption above 450 meV 
increases with current. As the modal gain g m is related 
to the material gain g via g m — Tg — a tot , where T is 
the confinement factor and a t ot the total losses averaged 
over the whole device [H[ , the experimental data and the 
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FIG. 1: Calculated energetically resolved electron concentra- 
tion and modulus square of 19 Wannier-Stark states for the 
central period. The upper and lower panels correspond, re- 
spectively, to low and high biases. 



calculated results should have the same spectral shape, 
but do not need to agree fully quantitatively. 

In Fig. [3] the calculated and measured current-voltage 
relation is shown. Although the slope on a logarithmic 
current-scale is similar there is still a substantial differ- 
ence between experiments and theory; almost one order 
of magnitude in current or 20 % in voltage. 



A. Gain Spectrum 

The reason for the bias dependence of the absorption 
peaks (see Fig. [5]) can be understood by an analysis of 
Fig. [1] For low bias, (top panel of Fig. [T]) , the carriers 
essentially occupy the ground state in the active region 
^ z < 20 nm, and the absorption peaks at 200, 280, 
370, and 420 meV, can be related to the transitions to 
the excited states in the active region. However, the ab- 
sorption peak at 370 meV is absent at zero current in our 
simulation. Here we have an almost vanishing dipole ma- 
trix element due to the approximate even parity of both 
the occupied state and the state at 430 meV in the up- 
per Fig.[T] As this peak is present in the experiment, the 
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FIG. 2: Top: Experimental modal gain estimated by trans- 
mission data as g = j- In T for TM polarisation. Gradually 
increased losses for the energies below about 160 meV indicate 
the long wavelength transmission edge for the laser waveguide. 
Bottom: Calculated material gain g. In experiments the tem- 
perature was 10 K for zero and low current and approximately 
80 K for the highest current. 
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FIG. 3: Calculated and measured current-voltage character- 
istic. While the voltage in the simulation is defined as the 
voltage drop per period times the number of periods, the ex- 
perimental values include the voltage drop over the contact 
regions as well. For small biases the calculated currents are 
not sufficiently precise for a reasonable display on a logarith- 
mic scale. 

real structure must differ somewhat. In addition there 
is some high energy- absorption by transition towards the 
high energy states above the band edges of the barriers. 
However the matrix elements between the well-localised 
ground state and these fairly delocalised states are small. 

For high bias and current, bottom of Fig. [2 the charge 
is transferred to the injector region. Thus the main ab- 
sorption peaks at 200, and 280 meV are bleached. The 
higher energy absorption peaks at 370, and 420 meV ex- 
hibit a non-monotonic behaviour with current, as the cor- 
responding excited levels mix with continuum states. 



The movement of charges to the injector leads to a 
significant increase of high energy absorption, since the 
states in the injector region are more extended and thus 
yield larger dipole moments for transitions to the contin- 
uum states. 



B. Spatial resolution of gain 

The location of charge density as a function of bias is 
explicitly visualised by the spatial resolved gain plotted 
in Fig. [4) which is calculated with the following expres- 
sion (see appendix [A"| . 

/ \ fe 2 \ - ip* a {z)d z il) {z) - d z ip* a {z)ij>p{z) 

2ce 0v /i7^ 2m e (z) 

x zp a {n a ~n p )— P ' a (3) 

(E a -Ep + hw) 1 + Tf j a /4 

Here d z ip(z) the derivative of ip(z) with respect to z. 

One can clearly see that most absorption/gain within 
the energy range 100 meV < E opt < 350 meV is located 
in the active region while the higher energy absorption 
follows the electron density position, i.e. localised in the 
active region for small currents and to the injector region 
for high currents. A single transition oscillates spatially 
between gain and absorption. However, integrating over 
z provides the total absorption from Fig. [2j The gain 
transition at 150 meV in the lower part Fig.|4]only shows 
gain contributions. This demonstrates the successful op- 
timisation of the dipole matrix clement of the lasing tran- 
sition during the design of the laser structure. 



C. High energy absorption 

As discussed above, the absorption for photon energies 
above E opt = 400 meV is a marker for the location of the 
electron charge, as the injector states are more extended 
compared to the ground state in the active region. Ac- 
cordingly the dipole transition matrix elements with the 
states from the continuum above the barriers are larger. 
While this trend is clear both in experiment and simula- 
tion, the detailed structure differs. 

We believe that this is due to two reasons. First, the 
calculated absorption in this range strongly depends on 
the number of states taken into account. This is shown in 
Fig. O where one can see, that the spectra converges with 
increasing number of states only for energies up to about 
500 meV. Second, the quantitative determination of these 
high energy states is less straightforward, as standard 
models such as effective mass or k • p models, are less 
accurate at these high energies. Other minima than the 
T-point might also come into play. 

To our knowledge, high energy absorption in QCLs has 
not been studied in detail. With the common approach 
to calculate the Wannier-Stark states directly for a QCL 
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FIG. 4: Spatial resolved gain for a low bias (top panel) and a 
high bias (bottom panel). Below 100 meV there is extremely 
strong absorption (on this scale) in the injector region. The 
conduction band profile in the plot is merely for spatial ori- 
entation, while its energy is arbitrary. 
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FIG. 5: Calculated spectrum for 1030 mA with different num- 
ber of Wannier States. The gain peak height is correctly de- 
scribed by the 14 states included in the transport simulation 
while the high energy absorption changed drastically when 
adding more high energy states. 



embedded in a finite box, the high energy states are diffi- 
cult to access, as their extended nature is not compatible 
with an artificial boundary. The difficulty is overcome by 
our approach starting with Wannier states [l^] , where the 
spatial cutoff is replaced by a cutoff in energy with re- 
spect to the local conduction band. While transitions to 
these states can shed light on the localisation of charges, 



we did not find any indication, that they are of relevance 
for transport properties or the gain spectrum around 150 
meV. Thus they can be safely neglected as long as only 
these standard properties are of interest. 



D. Discrepancies in current 

Comparing simulations and experiments we observe 
good overall agreement. However there is one main dis- 
crepancy: the calculated current is substantially larger 
than the experimental value. We attribute these results 
partly to the fact that the parameters used in calcula- 
tions slightly differ from the real sample parameters and 
partly to the limitations of the theoretical method used. 
The voltage drop over the device region and the effective 
electron concentration are the main less-known experi- 
mental parameters. 

The voltage drop Contact outside the QCL structure, 
such as contact regions and cladding layers may add up 
to about one volt or more, which could explain the differ- 
ence in the current-voltage characteristics, Fig. [3J where 



NFd + V ct 



while V s i 



NFd, where F 



is the average electric field over the device region and 
Nd = 1.83 /xm is the total length of the QCL structure. 
Also there is some uncertainty regarding the real dop- 
ing level that could be slightly higher or lower than the 
nominal value. 

Judging from Fig. [3] and assuming that many-body ef- 
fects are small (like electron-electron interaction and the 
Pauli principle) the relation 



(4) 



seems reasonable, where I is the current, n the electron 
concentration, and Fq is a constant. This suggest that a 
small error in the actual field can give a large error in the 
current, while an error in the electron concentration gives 
the same relative error in the current. We can therefore 
conclude that the discrepancy in current-voltage between 
experiments and simulation can be explained both by a 
small voltage drop over the contacts and cladding region 
and an uncertainty in the actual doping density. 

The experimental gain g m (w) in Fig. [2] is calculated 
from transmission data T(oj) with g(uj) — -i- lnT(co>), 
where L is the length of the laser ridge, 2.5 mm in this 
device. This value differs from the calculated material 
gain g due to the confinement factor and the frequency- 
dependent losses, which makes a quantitative comparison 
difficult. 

The same parameter uncertainty as in the current- 
voltage case affect to the overall magnitude of gain and 
absorption peaks, which are somewhat larger in the sim- 
ulations. If we assume, as previously, that many-body 
effects are small, a simple relation 



g(u) oc nf(F,u) 



(5) 



holds. Thus the onset of gain is related to a fixed field 
fonsct and for currents, where the gain spectrum exhibits 
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qualitative changes, quantitative information can be ex- 
tracted. In particular, the gain peak at 150 meV sets 
in at Texp ~ 600 mA, while this feature appears in the 
simulations for J s j m 1000 mA, which suggest a smaller 
effective electron density. Another aspect refers to the 
width of the gain/ absorption peaks, which are larger in 
the simulation. This indicates the presence of correla- 
tions in the scattering environment for the states involved 
[l5j . which we neglected in our simple model. 

A last remark concerns the position of the gain peak. 
The positions of the absorption peaks in simulation agree 
excellent with experiments, while in experiments the gain 
peak is 15 meV higher in energy than in simulation. This 
may be related to pronounced low frequency absorption 
features (such as free carrier absorption in the cladding 
regions) not taken into account in the simulation. Such a 
mechanism reduces the left side of the gain peak stronger 
than the right side thus shifting the observed gain peak 
position to higher frequencies. 



IV. CONCLUSION 



These features can be described by resolving the gain 
spatially, Eq. (J3]) . Our simulations are consistent with the 
experimental findings if we assume some uncertainty in 
the actual electron density in the sample and that a part 
of the voltage drops outside the actual QCL structure. 
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APPENDIX A: SPATIALLY RESOLVED GAIN 

In order to pin-point the spatial location of different 
spectroscopic features without resorting to investigating 
every possible transition an expression for a spatially re- 
solved gain is sought. Our starting point is the common 
relation fijl 



The evolution of the gain/ absorption spectrum of a 
quantum cascade laser provides a clear signal for the spa- 
tial transfer of charge inside the structure. If the elec- 
trons dominantly occupy a localised level (such as the 
ground state in the active region), one observes several 
distinct peaks due to transitions to higher levels in the 
same region, but only weak absorption to the delocalise 
states above the barriers. This continuum absorption be- 
comes more pronounced for transitions from the more ex- 
tended injector states, while the aforementioned absorp- 
tion peaks vanish, if the majority of charge is transfered 
hither. 



?R.{<t{uj,z)} 



ceo 



(Al) 



where we have taken into account a z-dependence of the 
gain g(uj,z) and conductivity cr(u>, z). The expression 
have a simple intuitive explanation: if the electric field 
and current are in phase, energy is accumulated in the 
sample corresponding to absorption (also known as Joule 
heating), and vice versa. 

The TM mode exhibits a small oscillating electric field, 
SF(lu), which is homogeneous on the length scale of a 
period of the structure. This implies the change 



for the Green's functions calculated in linear response from Ref. [l5| neglecting the ST, terms. The local current, 
5J(ui,z) is then given by [l3[ 



r 



Integrating over E and using a(ui,z) = 5Jq{lu,z)/8F(lu) 
we obtain Eq. ([3]). Averaging Eq. ([3]) over z we obtain 
Eq. (ffj by using 

^/^^^^W"^^^ (A4) 
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